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£T) • Abstract 

Ground-state properties of the two-dimensional 5 = 1/2 random Heisenberg models 
— |- . are investigated by the exact-diagonalization method. The phase diagram of the bond- 

random model (the ± J model) is the same as that of the corresponding Ising spin glass 
model, in spite of quantum fluctuation. In the site-random model, the spin glass phase 
exists at least in the ferro-rich case. In comparison with the corresponding Ising model, 
OS . the Neel phase of this model becomes much narrower in the ferro-rich case, while it is 

comparable in the antiferro-rich case. 
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While various aspects of quantum effects in random spin systems have been investigated, the 
recent progress in the random quantum Heisenberg model, the most realistic model for spin 
glass materials, is restricted to studies in one dimension. |], 0, |3| Quite recently, one of the 
present authors (Y. N.) investigated the 5 = 1/2 asymmetric Heisenberg Mattis model in 
two dimensions. |5|, He found that its ground-state phase diagram does not change^] from 
that of the corresponding Ising Mattis model, || in spite of quantum fluctuation. He also found 
that accurate calculations are possible only using quite small clusters (up to the 20-spin one) 
in the antiferromagnetic region. 

In classical spin glass systems, unexpected behavior was reported for the site-random Ising 
spin glass model in two dimensions || by Shirakura et a/.Q They analyzed this model using 
Monte Carlo simulations, and found that the spin glass phase may be stable at finite tempera- 
tures. The present authors investigated [H the ground state of the same model in detail using 



the numerical transfer- matrix method, |10j and found that the ferromagnetic and the Neel phase 
regions are wider than the expected values obtained from the corresponding ±J model. They 
also proposed that such structure of the phase diagram can be understood well in terms of the 
site-percolation picture. 
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On the basis of these findings, we investigate ground-state phase diagrams of several two- 
dimensional S — 1/2 random Heisenberg models described by the following Hamiltonian, 

H = JijSi ■ Sj . (1) 

(ij) 

Properties of the models depend on the choice of randomness in J^-. We consider the conven- 
tional bond-random model (the ± J model) [|TTJ in which the distribution function of is given 
by 

P(J tJ ; P ) = P 5(J i:j - 1) + (1 - p)S(J ij + 1), (2) 
and the site-random model 0, §, [| defined through 

Jij = -[e(l-w i w j )+Ui + u j ] , (3) 

where Ui is an independent random variable with the distribution 

P{ui) c) = c8{ui - 1) + (1 - c)6(ui + 1) . (4) 

The latter model is characterized by the Ising variable u>i located at each site. This parametriza- 
tion represents mixed crystals; i.e. u>i = +1 and uj^ = —1 correspond to different ions. In the 
case e = +1, = +J unless u>i = — 1 and uij = —1. We call this case "ferro-rich" (F- 
rich). Similarly, in the "antiferro-rich" (AF-rich) case (e = — 1), — — J unless — +1 and 
Dj = +1. Note that these two cases are not equivalent in quantum spin systems. 

We have two purposes in the present study. One is to clarify effects of quantum fluctuation in 
spin glass systems. Since quantum fluctuation plays an important role in the antiferromagnetic 
region, we calculate the Neel order parameter m st defined in 

m ^4EH) N, [(5r^] r . (5) 

Our other purpose is to investigate the difference between bond-random systems and site- 
random systems. As in the Ising spin glasses, we consider the spin-glass (SG) order parameter 
m sg defined in 



<4e[(^4) 2 L • (6) 

In Ising systems, the lower critical dimension d\ c of the SG phase has been expected to be 
in the range 2 < d\ c < 3[|l^] in the bond-random model, while it is expected to be equal to that 



of the pure Ising model, d\ c = l|] in the site-random model. In quantum Heisenberg systems, 



though d\ c has been reported to be greater than three fl3|, [L4J] in the bond-random model, the 
above argument suggests that it is equal to that in the pure quantum Heisenberg model, d\ c = 2, 
in the site-random model. If this is the case, the SG order of the three-dimensional site-random 
Heisenberg model is expected to be stable even at finite temperatures, which is consistent with 
experimental results. 



*If the existence of the SG order parameter at finite temperatures is established in two dimensions as 
pointed out by Shirakura et aZ.,@ d\ c — 1 is trivial. The one-dimensional nearest-neighbor spin glass model is 
not frustrated, and can be identified with the ferromagnetic Ising model by means of a gauge transformation. 
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We use the exact-diagonalization method in the present letter. Although the size of treatable 
clusters is quite small in this method, plausible estimation is expected in the antiferromagnetic 
region, as in the case of the asymmetric quantum Heisenberg Mattis model. j|] Calculations 
are carried out on the 8,10,16,18,20-spin Oitmaa-Betts-type||15|| clusters. In the site-random 
model, we average all the configurations. In the bond-random model, we average all the bond 
configurations in the 8 and 10-spin clusters, and average one thousand randomly chosen samples 
in the 16, 18 and 20-spin clusters. 

First, we show the data of the Neel order parameters. These quantities in the bond-random, 
the AF-rich site-random and the F-rich site-random models are plotted in Figs. |l|(a), |l](b) and 
|l](c), respectively. They are scaled well by iV -1 / 2 , as in the two-dimensional pure antifer- 
romagnetic Heisenberg model (p = or c = case).|16|] By the least-squares fitting of the 
10, 16, 18, 20-spin data, we estimate the critical concentrations of the Neel phase as 

p c = 0.11 ± 0.01 (bond-random model) , (7) 
c c = 0.36 ± 0.03 (AF-rich site-random model) , (8) 
c c = 0.13 ± 0.02 (F-rich site-random model) . (9) 

In the corresponding Ising spin glass models, these critical concentrations were estimated as 



Pc = 0.11 ± 0.01, [0, |c c = 0.37 ± 0.01 and c c = 0.41 ± 0.01,0 respectively. 

In the bond-random Ising model (±J Ising model), the Neel phase boundary in the p-T 
plane is vertical to the p-axis.JH], |18[ Namely, the Neel order is not destroyed by thermal 
fluctuation in the low-temperature region. If the effect of quantum fluctuation is similar to 
that of thermal fluctuation, the equivalence of p c in the Heisenberg model and the Ising model 
can be understood. Numerically, the estimates of c c are also consistent with each other in 
the AF-rich site-random Heisenberg model and Ising model. Although verticality of the phase 
boundary has been shown only in the bond-random model, a similar mechanism might work 
for the Neel order of the AF-rich site-random model. On the other hand, the estimates of c c 
are quite different in the F-rich site-random model (c c ~ 0.13 in the Heisenberg model, while 
c c ~ 0.41 in the Ising model). This nontrivial quantum effect cannot be explained by the above 
thermal-quantum equivalence. A study of the origin of such large difference is now in progress. 

Next, we show the data of the SG order parameters. These quantities are scaled by 1/N, 
as pointed out in the case of the two-dimensional Heisenberg Mattis model. @] An example of 
the fitting is given in Fig. |2| for the F-rich site-random model. Similar good scaling behavior is 
also observed in the other models. The estimates of the SG order parameter are plotted versus 
p or c in Fig. |[ These results suggest that the SG order parameter in the site-random model 
is nonvanishing in the whole parameter region of the concentration c in both the F-rich model 
and the AF-rich model . However, even in the bond-random model, the SG order parameter 
also seems to be nonvanishing. In the two-dimensional bond-random Ising model, this order 



parameter is stable only in the vicinity of p = and p = 1|]ID, |TT| in the ground state. This 
should also be the case in the present bond-random Heisenberg model, in which the SG order 
is weakened more by quantum fluctuation. Thus, there still remains a finite-size effect in the 
present calculation, and the existence of the SG order cannot be concluded only from these 
results. In fact, in the AF-rich site-random model, the estimate of the SG order parameter 
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decreases as c increases, and becomes as small as those in the bond-random model around 
c~ 0.5. 

The estimate of the SG order parameter in the Heisenberg Mattis model[|J is also displayed 
in Fig. |3|. The previous study of this modelQ] showed that the SG order parameter is nonva- 
nishing at all concentrations. The concentration dependence in the F-rich site-random model 
is similar to that in the Mattis model, except for small reduction around c ~ 0.5. Even in this 
region, the estimate of the order parameter is larger than that at c = 0. Thus, we can at least 
claim that the SG order remains nonvanishing in the F-rich site-random Heisenberg model in 
two dimensions. 

In summary, we study ground-state phase diagrams of various S — 1/2 Heisenberg spin 
glass models in two dimensions. In the bond-random model (the ±J model), the critical 
concentration of the Neel phase is estimated as p c = 0.11 ± 0.01, which is consistent with the 
corresponding Ising spin glass model. Such equivalence is expected to be related to the vertical 
phase boundary in the phase diagram of the bond-random Ising model. In the site-random 
model, the critical concentration is consistent with the Ising spin glass model in the AF-rich 
case (c c = 0.36 ± 0.03), while that in the F-rich case (c c = 0.13 ± 0.02) is decreased much 
more than that of the corresponding Ising model (c c ~ 0.41). This reduction is the only clear 
quantum effect displayed in the present study. Moreover, the evidence for the existence of the 
spin glass phase is found in the two-dimensional quantum Heisenberg site-random model in the 
ground state, at least in the F-rich case. Studies on the mechanism of the large reduction of the 
Neel phase in the F-rich site-random model and structures of excitation spectra in the present 
models are now in progress. 
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Figure Captions 



Figure 1: The Neel order parameters of 8,10,16,18,20-spin clusters are plotted versus N~ l l 2 
(a) in the bond-random Heisenberg model for p = 0.0, 0.1, 0.12, 0.15 (from above to below), (b) 
in the AF-rich site-random Heisenberg model for c = 0.0, 0.2, 0.3, 0.35, 0.4 and (c) in the F-rich 
site-random Heisenberg model for c = 0.0,0.05,0.1,0.15,0.2. The straight lines are drawn by 
the least-squares fitting of the last four data of each set. 

Figure 2: The SG order parameters of 8,10,16,18,20-spin clusters are plotted versus 1/N in 
the F-rich site-random Heisenberg model for c = 0.0,0.2,0.4,0.6 (from above to below). The 
straight lines are drawn by the least-squares fitting of the last four data. 

Figure 3: The estimates of the SG order parameter in the thermodynamic limit in the bond- 
random model (solid line with crosses), the AF-rich site-random model (solid line with circles), 
the F-rich site-random model (solid line with diamonds) and the asymmetric Mattis model 
(broken line). 
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